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Abstract 

The set of diagonals of real symmetric matrices of given non negative spec¬ 
trum is endowed with a measure which is obtained by the push forward of 
the Haar measure of the real orthogonal group. 

We prove that the Radon Nicodym derivation of this measure with respect 
to the relative Euclidean measure is approximated by the coefficients of a 
sequence of zonal sphere polynomials corresponding with the given spec¬ 
trum. There is a striking similarity between the role of the zonal sphere 
polynomials in the orthogonal case, and that of the Schur function in the 
Hermitian case. 

Following this we obtain a combinatorial approximation for the probability 
of real symmetric matrix of a given spectrum to appear as the sum of two 
real symmetric matrices, each of a given spectrum. In addition we obtain a 
real orthogonal analogue to the Zuber Itzykson Harish Chandra integration 
formula. 

Preliminaries and main result 

We use Greek letters to denote partitions of unity, 

6 = 5i,52---6n where di > 0 and partitions of 

integers, 7 h A. 

Let A = (Ai, A 2 ,...., A„) be a positive real vector so that ^ A* = 1 

For simplicity, we assume that all Xi are rational. 

Define a diagonal matrix by 

Dx = diag{Xi,X 2 ,...,Xn) (1) 
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Let On be the real orthogonal group of n by n matrices. 

By the Horn-Schur theorem, the set 

{{oDxO^)ll, {oDxO%2, •••, {oDxO^)nn/o G On) (2) 

is the permutahedron PHx, defined by the vector A. That is the convex hull 
of all the vectors (Ao-(i), Ao-( 2 )---Ao-(n)) for a G Sn- 

PHx is endowed with a measure, given by pushing forward of the Haar mea¬ 
sure of the real orthogonal group. We denote this measure by DH^, (In the 
Hermitian case DH^{see [FG])) 

We asumme that the Radon Nicodym derivation of DH^ with respect to 
the relative Euclidian measure of PHx exists, and is almost everywher con- 
tiuous. Furthermor the suspected points of the Radom Nicodym derivation 
function are included in severel afine spaces, their intersections with PHx, 
is of posetiv codimension. The proof will be given elswhere. 

Remark As not as in our case, the Radon Nicodym derivation of DH^ 
with respect to the relative Euclidian measure of PHx can be computed di¬ 
rectly using representation theory, thereby avoiding any need for differential 
geometry. See[F G] for more details. 

Let A be a partition of a positive integer N (A h N). Eurther,let 


Z,(X) = a nX^ (3) 

rjhN 

denote the zonal sphere polynomial corresponding to the partition A (see 
[Mac] section 7 and [J Ij) 

The next thorem is our main result. 

Theorem 1. Let A = Ai < A 2 ... < A^ he a non negative rational vector so 
that = 1, and let N be a positive integer, such that NX is a vector of 
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integers. 

Let Z]^x{X) = J2ri\-N corresponding zonal sphere polynomial. 

Assume ^ is a continuous point of the density function (Radon Nicodym 
derivation of DHf with respact to the relativ Euclidian measure on PH\). 
Than as N tends to infinity, corresponding unit matrixe) 

approximates this density function at the point ^ in PH\. 

Remark Compare to the characterization of the zonal spher polynomials 
coefficient in [Mac 2.27]page 409 

Remark In the Hermitian case the zonal sphere polynomials are re¬ 
placed by the Schur functions in order to obtain a similar approximation. 
See[F G] for more details 


Corollaries of the main result 

Theorem 1 leads to two corollaries we introduce now 

Let 7 = ( 71 , 72 , •••,7n) be another vector of rational non negative numbers 

such that — — diag{'yi,'y2, ■■■,'yn)- 

Let N be a large integer so that NX and N'j are integral vectors . Let 

ZNx{X)ZN^iX) = CY^^Z^iX) (4) 

ri\-2N 

be the linearization of the product of the two zonal sphere polynomials.(see 
[Mac] page 409 ) 


Claim 1. Using the Haar measure of the orthogonal group as a probability 
measure, the Radon-Nikodym derivative at ^ of the distribution of the ran¬ 
dom variable,^ satisfying the equation = D^/ ■. 0 , 0 ' ^ On, with 


N 

respect to the relative Euclidian measure on the simplex, is approximated, 
up to normalization, by _ 
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Simply speaking, the probability to obtain a matrix with spectra close to 
^;r] N, as the sum of two symmetric matrices of spectra A and 7 , chosen 
randomly with respect to the Haar measure of the real orthogonal group, 
can be evaluated by the linearization coefficients as above 

The reason for this is that one can define a measure on the set 

{liiD- + D;')n, {D'i + D°')22, {Dl + D-^U)} 

(using the notation A° = oAo^ ), given by the push forward the Haar mea¬ 
sure of On X On- Since this measure is DH^ * DH° (the convolution of the 
two measures),the main result tells us that it can be approximated by the 
coefficients of the product Zj^\{X)Zp^,y{X) 

Remark In the Hermitian case, the Littlewood Richardson coefficients re¬ 
place the coefficients in (4). See [F G] for more details 

The second corollary deals with the real symmetric version of the Zuber 
Itzykson Harish Chandra integral. 

Assume that C is a real symmetric matrix. For A and Dx as in ( 1 ), define 
IxiC) = j exp{trace{CD°x))dO (5) 

integrated over On with respect to its Haar measure. 

Without loss of generality, one can assume that C is a diagonal matrix. 
Now, the integral may be taken over PHx, with respect to DHy 
Ci,C 2 ---Cn be the spectra of C, then for arbitrary N we have 

hiC) = f (6) 

JPHx 

By using Riemann-sums to approximate the integration over PHx, comput¬ 
ing the integrand at the points of type ^ : 7 h A" in PHx, and then by 
the main result approximate the density function of DH^ at those points 
using the coefficients of the zonal sphere polynomial Z]xfx{X) we obtain the 
following. 
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Claim 2. If X is a unit partition as given above, then 

T tn\ ^ 7 ( Cl C 2 Cn -^ 

Remark One can ignore the points of PH\ in which the Radon Nicodym 
derivation is not continues thanks to our assumptions on DH\. 

Remark In the Hermitian version of this process, one may replace the zonal 
sphere polynomial in the right hand side with the Schur function with re¬ 
spect to the same parametrization. The very elegant determinant formula of 
Harish Chandra Zuber Itzykson is accepted as the limit of the determinants 
formula for the corresponding sequence of Schur functions. Unfortunately, 
the zonal sphere polynomial, do not yield such nice expression for their com¬ 
putations 


Proof of the main result 

Proof. Let us define a sequence of polynomials as follows: 

pf (X) = [ trace^{DlDx)do (8) 

J On 

once integrated over the orthogonal group with respect to its Haar measure, 
and X is a vector of indeterminants Xi,X 2 , ■■■Xn. 

Remark In [J 2] Identity 

Pf (X) = c,Z,{Dx)Z,{Dx) (9) 

r]\-N 

was obtained, and the coefficients Crj were investigated .However we avoid it 
in our treatment . 

Now we rewrite (8)as 

traced {DlDx) = (^{DDuXif (10) 

One can observe the similarity between the sequence of polynomials, P^ (X), 
and the sequence of Bernstein polynomials approximating the density of 
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DH^ (see [Lo] page 54) 

More explicitly, by writing 

= ( 11 ) 

r;l-Ar 

and combining with (10), we obtain 

<=( 12 ) 

i 

Where the integral is over the orthogonal group with respect to its Haar 
measure. 

Let Y = Yi, Y 2 , ...Yn be a real vector such that = ^ ■ Y 0. 

The next polynomials 

i 

are special members of the family of ” Dirichlet distributions”, so the inte¬ 
grand of equation (12) tends to 5^ (the delta function concentrated around 
the points over the n-1 dimensional simplex, as N tends to infinity 

Therefore, as N tends to infinity, one finds that the coefficients tend 
to the density of DH‘^ at the points ^ in which it is continuous. 

Next, we treat the left hand side of (10) by writing 

trace^{DlDx) = trace{{DlDxf^) (14) 

Let V denotes an n dimensional vector space over the reals, and let vi,V 2 , --Vn 
be an orthogonal basis of V. 

Let us define the r] weight space for any partition r/ of N, of no more than n 
parts, by the formula 

= span{vi^ (g) Vi2 0 ... 0 Vij^ : #{/ ; ii = j} = r]j} (15) 
Now because Dx is a diagonal matrix, we have 
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trace^ {DxDx) = trace{{DxDx)®^) = '^trace{D\Dx)®^)\v'SN^ (16) 

V 


and 


traceUDxDxnv^N^ " ( ^ J 


(17) 


Now, using spectral factorization of real symmetric operators, one obtain 
that 

( 0 ,)®" = ( 18 ) 

?7l-Af 

where the second summation runs over a basis of the the eigenspace of weight 
7] with eigenvalue 7r? = 0 ) 

(One can choose the basis as the standard unit vectors in V^^). 


Now,one can use the central limit law (using equation (17)) to write 
Px{X)r^ f jjj y^ trace{{o‘^^a^){o‘^^a.yYD^^)do (19) 

On AT 


r]\-N 


summed over rj}- N such that ^ ~ A. 

The first ~ in the equation means that contribution of the other summands 
to the coefficients vanishes to zero. The second ~ denotes the closeness 
of real vectors at any coordinate. 

To get closer to our goal one may apply representation theory of Gln{R) 

By the Schur Weyl duality we have 


= 0 ® M,, (20) 

where the sum is over rjh N of no more than n parts, 

and Vrf{Mrj) are the GlniR){Sx) simple modules in their actions on the 
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tensor product space.See [G W] page 138 for more details 
The striking fact is that, as N tends to infinity 

trace^{DxDx) ~ ^ trace{{DxDx)^^\vr,^Mr,) (21) 

(see the remarks on~ after (19)). 

The discovery of this concentration phenomena is, some times, attributed 
to Kyel and Werner and has been rediscovered several times. (See [Ch] for 
a nice proof of it). 

Remark At this point, one can easily see the full Hermitian perspective of 
the problem. See [F G] for more details . For symmetric real matrices one 
has more work to do. 

Using (19) and (21), one can write 

Pf(X)~ f Y,c^itrace{{o^^a^){o^^a^YDf^))do (22) 

7 

summed over a basis of eigenvectors (of eigenvalues Cy) of in (U 5 / (8> 

My)n(^®%,where^,|;~A. 

Now, for each vector a in , we define a polynomial by the formula 
Qa{m) = trace{a^am®^) (23) 


for any n by n matrix m 

Before introducing the next lemma, we quote the definition of zonal sphere 
polynomials from the first page of [J 2] 

”... Denote by Vj the space of polynomials of degree f over the real sym¬ 
metric positive definite n by n matrices. 

GLn{R) acts onVf as follows 

(For g e GLn{R), P ^Vf : giPim)) = P{gmg*)). 

As GLn{R) module Vf = 0(j) F(j) where (/) = fi, f 2 ---fn is a partition of 



f and is the subspace of Vf on which GLn{R) acts irreducibly corre¬ 
sponding to the partition 2(f) 

The zonal polynomial Zj spans the real orthogonal group invariant vector 
space in Vj... ” 

To complete the proof of Theorem 1 we need the following lemma 

Lemma 1. Using the notations from [J2], for as in (22),we have that 
Qa^im) G 0b);/) summed over (f) such that ^ ~ A. 

Proof. One need to apply some basic representation theory of gln{R), the 
Lie algebra of Gln{R) ■ [G W] page 228 is good reference to refer. We note 
that the weights referred to in (15) and the weights in the context of Lie 
algebras representation theory coincide for an appropriates choice of basic 
positive roots vectors. 

Now, since the map a i—)■ of a is bilinear, it can be factorized through 
the tensor product, a i—)• a 0 a i—)• of a. 

Furthermore ior 5 h N we haveV^ 0 = 0 Kj, summed over lo ^ 25, 

where ^ is the dominant order defined on the weight lattice ( see page 237 
in [G W]) 

Since a.y is of weight 5f 0-^0 a.y is of weight 25', and so belongs to 0 
summed on uj in the next small fragment; 2N(X — e) ^ 25' :< tv P 2N{X + €'), 
where e, e' are arbitrarily small positive vector , as N tends to infinity . The 
central inequality just expressed the fact that the weights of vectors in a 
highest weight modules are smaller than its highest weight. □ 

Given Lemma 1, Theorem 1 is now proved since the integral in (22) is 
just a projection onto the space spanned by the unique zonal polynomial in 
each V)/). 

□ 
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Numerical Experiments: 


We test the 3x3 case for the diagonal matrix dza(7(0,1,2), conjugated 
by (a) SO{3), and (b) SU{3) (Figure 1). In each case we sample 30,000 
conjugations, and compute the resulting diagonal. In both cases we have 
a two-dimensional distribution, based on the (1,1) and (2,2) positions of 
the diagonal. The (3, 3) entry is determined from the first two by the trace 
condition. . 

We have also test conjugations of rank 1 matrices 

We sample the diagonals of a 3 x 3 matrix of rank 1 conjugated by SO{3) 
random matrices (See Fig 2 (a)) On the other hand we sample the zonal 
sphere polynomial’s coefficients corresponding to one row (using example 1 
of[Mac]Page 410)(See Fig 2 (b)) 

In Figure 3 the coefficients of the zonal sphere polynomial of two variables 
corresponding to one row (the red points), have been compared with the 
function; xy~ 2 , x + y = constant, 

which is computed in [FG] explicitly for two by two symmetric matrices for 
that case. In the figures, we use the notation as employed in [Mac, p. 
385] for the Jack-polynomials. Note that the r = 2 case corresponds to the 
zonal polynomials. 



(a) Conjugation by SO{3) (b) Conjugation by 5(7(3) 

Figure 1: 
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Figure 2: 



Figure 3: Zonal coefficients of J|q (dotted) against theoretical prediction 
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